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Abstract 

We define and study on Lorentz manifolds a family of covariant diffusions in 
which the quadratic variation is locally determined by the curvature. This allows the 
interpretation of the diffusion effect on a particle by its interaction with the ambient 
space-time. We will focus on the case of warped products, especially Robertson- 
Walker manifolds, and analyse their asymptotic behaviour in the case of Einstein-de 



Sitter-like manifolds. 
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1 Introduction 



It is known since Dudley's pioneer's work [Du] that a relativistic diffusion, i.e. a Lorentz- 
covariant Markov diffusion process, cannot exist on the base space, even in the Minkowski 
framework of special relativity, but possibly makes sense at the level of the tangent bundle. 
In this spirit, the general case of a Lorentz manifold was first investigated in [F-LJ], where 
a general relativistic diffusion was introduced. The quadratic variation of this diffusion is 
constant and does not vanish in the vacuum. In this article, we investigate diffusions whose 
quadratic variation is locally determined by the curvature of the space, and then vanishes 
in empty (or at least flat) regions. 

The relativistic diffusion considered in [F-LJ] lives on the pseudo-unit tangent bundle 
T X M of the given generic Lorentz manifold (A4,g). As is recalled in Section 3.1 below, it 
can be obtained by superposing, to the geodesic flow of T Y M., random fluctuations of the 
velocity that are given by hyperbolic Brownian motion, if one identifies the tangent space 
T^Ai with the hyperbolic space M d (at point £ G M., by means of the pseudo-metric g). 
These Brownian fluctuations of the velocity can be defined by the vertical Dirichlet form 

/ V^F(£,£) n(d£,d£) , considered with respect to the Liouville measure fi . 

The Dirichlet forms we investigate in this article depend only on the local geometry of 
(A4,g), e.g. on the curvature tensor at the current point £, and on the velocity £. We 
consider several examples : 

- If the scalar curvature R(£) is everywhere non-positive (which is physically relevant, see 
[L-L]), then the Dirichlet form can be : / V"F{£, £) 2 R(£) n(d£, d£) , 

JT X M 

leading to the covariant relativistic diffusion we call R- diffusion . 

- If the energy £(£,£) is everywhere non-negative (which is physically relevant, see [L- 
L], and [H-E], where this is called the "weak energy condition"), then we can choose the 

Dirichlet form to be : / V£F(£, £) * S(£, £) d£) , 

leading to another covariant relativistic diffusion, we call energy ^-diffusion . 

Contrary to the basic relativistic diffusion, these new relativistic diffusions reduce to the 
geodesic flow in every empty (vacuum) region. 

Note that —R(£) and £(£,£) could be replaced by (p(R(£)) and ip(£ (£,£), for more 
or less arbitrary non-negative increasing functions tp, if). We shall present this class of 
covariant S-relativistic diffusions , or S-diffusions , in Section 3 below. 
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- If the sectional curvatures of timelike planes are everywhere non-negative (sectional cur- 
vature has proved to be a natural tool in Lorentzian geometry, see for example [H], [H-R]), 
as this is often the case (at least in usual symmetrical examples), then it is possible to 
construct a covariant sectional relativistic diffusion which undergoes velocity fluctuations 
that are no longer isotropically Brownian, using the whole curvature tensor (not the Ricci 
tensor alone). See Section 6 below. This sectional relativistic diffusion depends on the 
curvature tensor in a canonical way. Its diffusion symbol vanishes in flat regions (i.e. re- 
gions where the whole curvature tensor vanishes), but does not vanish, in general, in empty 
regions (i.e. regions where the Ricci tensor vanishes). 

Note that all these covariant diffusions are the projections on T X M. of diffusions on the 
frame bundle G(M). Actually, they are constructed directly on G(Ai), as in [F-LJ] and 
in the classical construction of Brownian motion on Riemannian manifolds, see [El], [M], 
[I-W], [Em], [Hs], [A-C-T]. These constructions are performed in Sections 3 and 6 below. 

Note also that, while in the flat case the Dudley diffusion [Du] is the unique covariant 
diffusion, the above examples show that this is not at all the same for curved spaces. 

In Sections 4 and 5, we study in more detail the case of warped products, and spec- 
ify further some particularly symmetrical examples, namely Robertson- Walker manifolds, 
which are warped products with energy-momentum tensor of perfect fluid type. 

We investigate more closely Einstein-de Sitter-like manifolds (Robertson- Walker mani- 
folds for which the expansion rate is a(t) = t c for some positive c), reviewing in this simple 
class of examples, the relativistic diffusions we introduced, which appear to be distinct. 
We perform in this setting an asymptotic study of the /^-diffusion, and of the minimal 
sub-diffusion (t s ,t s ) relative to the energy ^-diffusion. 



2 Canonical vector fields and curvature 

We present in this section the main notations and recall a few known facts (see [K-N]). 
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2.1 Isomorphism between /\ E ' and so(l,d) 



On the Minkowski space-time IR 1,d , we denote by rj = ([rfijl)o<i,j<d the Minkowski tensor 
/l ... \ 
-1 ... 

. We also denote by ((?7 y ))o<i,j<(i the inverse tensor, so that JfoV = 5f 

V° ••• -V 

(or equivalently : 77^ = rfi := l{j=j=o} — l{i<i=?<d}), an d by (•,•),, the corresponding 
Minkowski pseudo-metric. For u,v,w G lR 1,a! , we set 



uAv(w) := (u,w) v v—(v,w) v u. 



(1) 
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In other terms, this is the interior product of u A v by the dual of w with respect to i] . 

This defines an endomorphism of R l,d which belongs to so(l, d), since for any w, w' G M 1 ' d 
we have clearly (u A v(w),w') r) + (w,u A v(w')) r] = 0. It vanishes only if u and v are 
collinear, hence if and only if u A v = . We have thus an isomorphism between f\ 2 M}' d 
and so(l,d). 

Remark 2.1.1 The Lie bracket of so(l,d) can be expressed, for any a,b,u,v G M}' d , by: 
[a A b, u A v] = (a, u) v b A v + (b, v) v a Au — (a, v) v b A m — (b, u) v a Av . 

The Minkowski pseudo-metric (•, -) v extends to /\ 2 IR 1 ' a! , by setting: 

(u A v, a A b) v := {u, a) v {v, b) v - (u, b) v (v, a) v = \{{u A v(a),b) v - (u A v(b), a) v ), (2) 

so that, if (eo, • • • , e^) is a Lorentz (i.e. pseudo-orthonormal) basis of (IR 1 '^,^), then 
{ci A tj | < i < j < d) is an orthogonal basis of (/\ 2 M}' d , r/), such that (e^ A €j , A e^)^ = 
Vu Vjj ■ 

2.2 Frame bundle G(M) over (A^,p) 

Let Ai be a time-oriented C°° (l+(i)-dimensional Lorentz manifold, with pseudo-metric 
g having signature (+, — ), and let T l M. denote the positive half of the pseudo- 

unit tangent bundle. Let G(A4) be the bundle of direct pseudo-orthonormal frames, with 
first element in T x hA , which has its fibers modelled on the special Lorentz group. Let 
i\\ \ u i — y (tt(m), eo(w)) denote the canonical projection from G(Ai) onto the unit tangent 
bundle T l M , which to each frame (e (w), . . . , e d (u)) associates its first vector e . 

We denote by TM^M the tangent bundle, by Y{TM) the set of C 2 vector fields on 
M (sections of 7T 2 ), by G(M) M the frame bundle, by u = (n(u); e (u), . . . , ed{u)) 
the generic element of G(Ai). We extend (1) to a linear action of so(l,d) = /\ 2 'R 1 ' d on 
G(M), by setting: 

e k A e £ (ej(u)) := r] jk e e (u) - r) je e k (u) , for any < j, k, I < d , 

where (e , . . . , e d ) denotes the canonical basis of M 1,d . 

The action of SO(d) on (d, . . . , e<f) induces the identification T X M = G(M)/ SO(d). 

The right action of so(l,d) on G(M) defines a linear map V from so(l,d) into vector 
fields on G(M) (i.e. sections of the canonical projection of TG(M) on G{M)), such that 

[V aA b, V aA p] = V [aAb>aAl3] , for any a A b , a A (3 G /\ 2 (3) 
Vector fields V aA b are called vertical . 

Notation To abreviate the notations, we shall consider mostly the canonical vector fields : 

Vij := V e ,Ae 3 , for < i, j < d . 
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By (3) and (1), for < i, j, k, £ < d we have : 

[% V ke ] = rj ik Vjt + rj je V ik - rj ie V jk - r] jk V u . (4) 

We shall often write Vj for V 0j . 

Denote by p the canonical projection TTM — TM , by p the canonical projection 
TG(M) G(M), and consider also the projection TTM {p ' T * 2 \ TA"! © , where 

TM®TM := {(£;v 1 ,v 2 )\€ e M ,v t e T t M } = { (w 1: w 2 ) \ Wi eTM, 7r 2 ( Wl ) = ir 2 (w 2 ) } 

is the so-called Whitney sum. 

A connection a can be defined as a bilinear section TM@TM — > TTM of (p,Tir 2 ), 
the bilinearity being that of (^1,^2) >->■ cr(C ',vi,v 2 ), above any given base point £ e 

Given such a connection cr and a C 1 curve 7 , the parallel transport //J along 7 of 
any t> G T 7o .M is t> t = //Jv E T Jt A4 defined by the ordinary differential equation : 
j|(7t ; v t ) = cr(7t ; Ut,7t). Then the covariant derivative V-y X(7 ) G T 10 M of a C 1 
vector field X is defined as the derivative at of t \-t (//J) 1 X('y t ) . 

A connection a is said to be metric if the associated parallel transport preserves the 
pseudo- metric, and then acts on G(M) as well as on TM. . 

A metric connection a defines the horizontal vector fields H k on G(M), for < k < d , 
given for any F G C 1 (G(M)) and u G G(M), by: 

H k F{u) is the derivative at of t \-+ F (fpu) , (5) 

the C 1 curve 7 being such that 70 = 7r(w),7o = e k {u) . Note that Tn(H k ) = e k . 

The canonical vectors Vij,H k span TG(M) (the horizontal (resp. vertical) sub-bundle 
of TG(M) being spanned by H k s (resp. V^'s)). Note that H generates the geodesic 
flow, that Vi, . . . , Vd generate the boosts , and that the Vij (1 < i,j < d) generate rotations. 

This allows to define the intrinsic torsion tensor ((7^)) and curvature tensor (iTZij ke )) (with 
< k,£ < d) of the metric connection a, by the assignment : 

d 

[Hi, Hj] =J2%*Hk+ E n ^ ^ ( 6 ) 

k=0 0<k<£<d 

we can denote more simply by l~j H k + \ 1Zij M V k e ■ 
For any metric connection we have : 

[V^, H k ] = 7] ik Hj - 7] jk Hi , for < i, j, k < d . (7) 
There exists a unique metric connection with vanishing torsion, called the 
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Levi-Civita connection . We shall henceforth consider this one. 
The curvature operator is defined on /\ 2 T^M. by : 

TZ^(ei(u)Aej(u)) := ^ e fc (u) Ae^(u), for any u G 7T _1 (^) and 0<i,j<d. (8) 

0<£i<f<(i 

The curvature operator is alternatively given by : for any C 1 vector fields X, Y, Z,A, 

(K (X A Y) , A A Z) v = ( ([Vx , V y ] - V [X , Y] ) Z,A) g . (9) 
The Ricci tensor and Ricci operator are defined, for < i, k < d , by : 

d d 

itf := ^ Tlij kj , and Ricci e (e^)) := R* e k (u) , for any u G 7r _1 (^) . (10) 
i=o fc=o 

The scalar curvature is : R := ^ -R^ . 

fc=0 

The indexes of the curvature tensor ((T^ij )) and of the Ricci tensor )) are lowered or 
raised by means of the Minkowski tensor ([i] a b)) and its inverse ((f] ab )). For example, we 
have : TZ p jqr = rf v i] kg r) ir , and R tj = R* r) kj . 

Remark 2.2.1 The curvature and Ricci operators and tensors are symmetrical: 

(lZ(a A b), v A w) v — (a A b, 1Z(v A w)) v , and (Ricci(t> ), w) v = (v , Ricci(w)) r? , 
for any a, b,v,w G Equivalently, for < k,£ < d: Rij U = R M ij, and Rij = Rji . 

The energy-momentum tensor ((X*)) and operator T$ are defined as : 

Tj :=R)-\R 5* and T 5 := Riccig - \ R . (11) 

d 

Note that Tj = — R . The energy at any line-element (£, £) G T X A^ is 

i=o 

The last equality is easily derived from (10) and (11) since writing (£, £) = (ir(u), e (u)) 
for any it G 7rf 1 (^,0 and Ty = T t k r/ k j , we have : 

( T t(£)>t)g(0 = 9{T^(e ), e ) = #(T fc e fc , e ) = T fe % = T 00 = T 00 (£, . 

The weak energy condition (see [H-E]) stipulates that £(£,£) > on the whole T X M.. 

We shall need the following general computation rule. 
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Lemma 2.2.2 For < i, j, k, £,p, q < d , we have : 

V 1Z ke = n 1Z M -n- TZ ke + n 1Z M -n- 1Z k£ +S k lZ- l -8 k TZ- e -S e lZ- k + 8 i 1Z- k 

v qp'^lj <{q% /v p.J 'lip /V 9J '113 /V *P UP /v «7 ' q /V «JP u p ,x "ijq u q /v «jp ' u p' x 'tjq • 

Proof Using (8) and (1), we have indeed: 

Vqi ^pj Vip ^qj Vqj ^ip Vjp ^iq ^q ^ijp $p ^ijq ^q ^ijp ^p^ijq • ^ 



2.3 Expressions in local coordinates 

d 

Consider local coordinates (£ l , e k ) for u — (£, e , . . . , e<i) G G(M), with Cj = e k ^— ^ ■ 

Then the horizontal and vertical vector fields Vij,Vj, H k , which satify the commutation 
relations (4), (6), (7) of the preceding section 2.2, read as follows. First, denoting by T\.- = 
F e jk the Christoffel coefficients of the Levi-Civita connexion V, we have for < i, j < d: 



v — = r k (n — 
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k m-ni 

'3 C * kmV*) q l 



(13) 



This is consistent with (5). Indeed, we have a priori an expression Hj = a k -J^ + > 

with on one hand a k = (Tir(Hj), d£ k ) = (ej, d£ k ) = e k . On the other hand, for a C 1 curve 
7 satisfying 70 = £ , 70 = ej , denoting by e _1 the matrix inverse to e = ((ef)) we have : 



^j x km Qti 
d n 



dt I 



_ d ■ d 

^70 ~XCm~ = ^0 (( e )m 6 i) = ^ 



d 







i^U//»\^ = Hj^Tm) 4jjg = -(e-^Hjiel) ^ 



by (5), so that = Hj(ef) = —e k e™T km as wanted. 

Recall that the Christoffel coefficients of V are computed by : 

r * = 1 n u( d 9i3 , 9gu d 9lj \ 
13 2 y V d? d& ) ' 

or equivalently, by the fact that geodesies solve £ fe + T k -^ 1 ^ 3 = . 
Then for < i, j, k < d : 



Vij e k = r) ik ej - T) jk ei = (r) ik e] - r) jk e\ ) 



d_ 



(Viq e j Vjq e i ) 



d , d 



1 ! de™ k dC~ 



whence for < i,j < d: 



V i3 = (vnef-mieT) 



_d_ 



(14) 



The curvature operator expresses in a local chart as : for < m, n, p, q < d , 

F r r s _-nr r s , OL nq np \ /-, r\ 

Then, the Ricci operator can be computed similarly, as : for < m, p < d , 

/ \ ~ dT n BY n 

J^mp ■ \^^^\ d ^ m ) > Q£ P J ^ l^mnpqy L nq L mp L pq L mn^~ ' 

The scalar curvature and the energy-momentum operator can be computed by : 

R = Rij g lj and T fm = R hn — \Rgt m . (17) 

To summarize, the Riemann curvature tensor ((Jlij ke )) is made of the coordinates of 
the curvature operator 1Z in an orthonormal moving frame, and its indexes are lowered 
or raised by means of the Minkowski tensor (( r) a bj), while the curvature tensor (lR mnpq )) 
is made of the coordinates of the curvature operator in a local chart, and its indexes are 
lowered or raised by means of the metric tensor ([g a b))- 

To go from one tensor to the other, note that by (15) and (8) we have 

^(|rA|r) = ^ n mn ab £ A ^ , whence : e* e] TZ M ^ = TZ l3 mn e? m e£ , or equivalently : 

n ijab = TZ Mrs e\ e] e r a e s b , or as well : n rs ™ = TZ abmn e r a e s b e p m . (18) 

2.4 Case of a perfect fluid 

The energy-momentum tensor T (of (11), or equivalently T, recall (17)) is associated 
to a perfect fluid (see [H-E]) if it has the form : 

Tm = qU k Ue-pg k £, (19) 

for some C 1 field U in T l Ai (which represents the velocity of the fluid), and some C l 
functions p, q on M. . By Einstein Equations (17), (19) is equivalent to: 

Rki = qU k U e + pg ke , with p = (2p-q)/(d-l), (20) 

or as well, by (16), to : 

(fficti(W), W) v = qx g(U, W) 2 + p x g(W, W) , for any W G TM . (21) 

The quantity (£/(&,), £ s ), (which is the hyperbolic cosine of the distance, on the unit 
hyperboloid at £ s identified with the hyperbolic space between the space-time velocities 
of the fluid and the path) will be denoted by A s or A(£, s ,t; s ). 

Note that necessarily A s > 1 . By Formulas (12) and (19), the energy equals: 

e(U) = q(Z)A(Z,£) 2 -p(Z)- (22) 
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Remark 2.4.1 (i) The energy of the fluid is simply: T k eU k U e = q — p. and the 
scalar curvature equals R — 2 [(d + 1) p — q]/(d — 1). 

(ii) By (22), the weak energy condition reads here : q > p + . 

3 Covariant 2-relativistic diffusions 

Let S denote a non-negative smooth function on G(A4), invariant under the right 
action of SO(d) (so that it identifies with a function on T X A^). 

Our examples will be S = — g 2 R and S = g 2 £ (for a positive constant g). 

We call S-relativistic diffusion or S-diffusion the G(.M)-valued diffusion process ($ s ) 
associated to S we will construct in Section 3.2 below, as well as its T 1 A^-valued projection 
7Ti($ s ). Let us first recall our previous construction, which corresponds to a constant S . 

3.1 The basic relativistic diffusion 

The relativistic diffusion process (£ s ,£ s ) was defined in [F-LJ], as the projection un- 
der 7Ti of the G(.M)-valued diffusion s ) solving the following Stratonovitch stochastic 
differential equation (for a given Revalued Brownian motion (w 3 s ) and some fixed g > 0) : 

d 

d^ s = H (V S ) ds + gJ2Vj(^s) o dwi . (23) 

3=1 

The infinitesimal generator of the G(A^)-valued relativistic diffusion (ty s ) is 

-H:=H + £J2V 2 , (24) 

3=1 

and the infinitesimal generator of the relativistic diffusion (£ s , £ s ) := 7Ti(\I/ s ) is the relativis- 
tic operator : 

■= ^ + i ^ = ^ + (¥ i k - ?£ s rj(o) ^ + 4 - 9 ke (0) ^ > (25) 

where C° denotes the vector field on T l Ai generating the geodesic flow, and denotes 
the vertical Laplacian, i.e. the Laplacian on T^Ai equipped with the hyperbolic metric 
induced by g(£). The relativistic diffusion (£ s ,£ s ) is parametrized by proper time s > 0, 
possibly till some positive explosion time. 

In local coordinates e*), setting \1> S = (£*,e*(s)), Equation (23) becomes locally equiv- 
alent to the following system of ltd equations : 

d 

= is ds = e k (s) ds ; d^ = -F k e (Q g 4' ds + g^ e?(a) dw\ + d -f g ds , and 

i=l 
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de%s) = -I* (&) e £ j(s) £ds + Q 4 fc dw j s + £ ej(s) ds , for 1 < j < d , < k < d . 
Remark 3.1.1 We have on T X M : 

d 

^V}£ = 2(d+l)£-2Tr(T) = 2(d + 1) £ + (d - 1) R . (26) 

3=1 

Indeed, since for each 1 < j < d Vj exchanges the basis vectors e = £ and ej (recall 
( 14) ) we get : Vf£ = 2 V 3 (f tm ^ ef) = 2 f tm (^ £ m + ej ef) , whence 

d 

VfS = 2d£ + 2T £m (i e r - g im ) = 2(d + 1) £ — 2 Tr(T) , 
and Formulas (26) follow at once by (11). 

As an application, a direct computation yields the following evolution of the energy. 

Remark 3.1.2 The random energy process £ s = £(£ s ,£ s ) associated to the basic relativis- 
tic diffusion tti(^/ s ) = (£ s ,£ s ) satisfies the following equation (where V„ := v 3 Vj) : 

d£ s = V is £ s ds + q 2 [(d + l)£ s + ^ R{Q] ds + dM £ s , 

with the quadratic variation of its martingale part dM £ given by : 

[d£ s , d£ s ] = [dMf, dM £ s ] = V [£ 2 S - (ft, ft)} ds . 

(We have here in particular = [%^-fe) - 2 f u {Q Y%{Q ] 

Note that the energy £ s is not, in general, a Markov process. 

3.2 Construction of the S-diffusion 

Let us start with the following Stratonovitch stochastic differential equation on G{M) 
(for a given Revalued Brownian motion (w 3 s )) : 

d d 

rf$ s = H ($ s ) ds+\J2 V 3 E (**)VA*s) ds + Y, V^W*,) o dwi . (27) 

3=1 3=1 



Note that all coefficients in this equation are clearly smooth, except vS on its vanishing 
set S _1 (0). However, y/E is a locally Lipschitz function; see ([I-W], Proposition IV. 6. 2). 
Hence, Equation (27) does define a unique G(.M)-valued diffusion (3> s ). 

We have the following proposition, defining the S-relativistic diffusion (or S-diffusion ) 
($ s ) on G{M), and (£ s ,£ s ) on T X M., possibly till some positive explosion time. 
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Proposition 3.2.1 The Stratonovitch stochastic differential equation (27) has a unique 
solution (<& s ) = (£ s ; £ s , ei(s), . . . , e^(s)) ; possibly defined till some positive explosion time. 
This is a G(Ai) -valued covariant diffusion process, with generator 

d 

n s := fffl+ij]^. (28) 
j=i 

Its projection 7Ti($ s ) = (£ a ,£ s ) defines a covariant diffusion onT x M., with S 0(d) -invariant 
generator 

-Hi := £° + ±V v ~V v , (29) 
denoting the gradient on T^M. equipped with the hyperbolic metric induced by <?(£). 
Moreover, the adjoint of H s with respect to the Liouville measure of G(M) is 

d 

%~ := —H + | ^2 Vj S Vj . In particular, if there is no explosion, then the Liouville 

3=1 

measure is invariant. Furthermore, if S does not depend on i.e. is a function on Ai, 
then the Liouville measure is preserved by the stochastic flow defined by Equation (27). 

We specify at once how this looks in a local chart, before giving a proof for both statements. 

Corollary 3.2.2 The T l Ai-valued E-diffusion (£ a ,£ s ) satisfies d^ s = £ s ds , and in any 
local chart, the following ltd stochastic differential equations: for < k < d, (denoting 
S s = S(£ s ,g)j 

dg = dM k s - r* (&) g g ds + lE s g ds + \ [g g - g k %)} i) ds , (30) 

with the quadratic covariation matrix of the martingale term (dM s ) given by: 

[dg, dg] = [g 4' - 9 k %)\ Zsds, for 0<kj<d. 

Proof In local coordinates (£', e k ), $ = (£ ; e , e±, . . . , e^), using Section 2.3, Equation 
(27) reads : for any < fc < d , dg = g = e k (s) ds , 



dg = -r&(6) S g + i E v i E fa> &) e % s ) ds + J2 V & e % s ) ° d < ; 

and for 1 < j < d , < k < d , 



de*( a ) = -I* (&) eS(a) g ds + J ^ S g) g ds + ^H(6,g)g ° . 



We compute now the Ito corrections, which involve the partial derivatives of yS(£, £) with 
respect to £ . We get successively, for l<j<d,0<k<d: 
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d V£(6,6K(s) ),dwi 



■(6,6) [dej(s),d^]+i 



2 eJ(s) 



= 2 (6, 6) 3/2 C fc ds + i e }( s ) -g(6, a^(C.C) e]{s) ds , 

OK 



hence, 



d h/S(6,6)e^) ),dw> 



s(6,6)4 fe ^ + i^5(6,6)e?( S )^; 



s(6,6 



S(6,6) dg,dw>. 



~ 2 



dZ(£ s ,£ s ),dw J s 



= S (6, 6) 3/2 e}( a ) ds + | £ ^(6, CSyHiCl! efc) ds , 



hence, 



d(VH(6,6)g),^ 



s(6,6)ej( s )^ + |^s(6,6)4 fe ^- 



Note that the simplification by yS(6,6) * s allowed, since on S 1 (0) both sides of the 
simplified formula vanish identically. 

Hence, in local coordinates and in ltd form, Equation (27) reads : for any < k < d , 



d£s = 6* ds = eg(s) ds , and setting dM s fc := £ J~(6, 6) ej(s) du?^ , 



3 = 1 



di k s = dM k s - r* (6) g g d s + 1 s (6,6) 4 fe ds + \ [g g - <^(6)] ^(6, 6) <fe ; 



de*(a) = V5(6, 6) t* dwi - I* (&) e (s) f> + ± E (6, 6) e}(a) ds + § ^ S (6, 6) C <** • 



Note that we used the formula e j( s ) e j( s ) = £s£f ~~ Q M {is) , which expresses that 

$ G G(.M). Using again this formula, we get the quadratic covariation matrix of the 
martingale term (dM s ), displayed in the above proposition, which shows that (7Ti($ s )) is 
indeed a diffusion, and proves Corollary 3.2.2. 

On the other hand, comparing the above equations with Equations (23) and (24), which 
correspond to S = g 2 , we get precisely the wanted form (28) for the generator of ($ s ). 

Then, since S = S A u + ( V* 1 S) , comparing the above equation (30) for 
with Equation (25) (for which S = g 2 ), we see that establishing Formula (29) giving the 

dE d 

projected generator "Hi reduces now to proving that (V S) = 6 — g ke ] — 7- — — • 

(K ( d£ k 
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Now this becomes clear, noting that V" = e k , — — , for each 1 < j < d . 

& 3 3 Qtk 

Finally, the assertions relative to the Liouville measure are direct consequences of the 
fact that the vectors H and Vj are antisymmetric with respect to the Liouville measure of 
G(A4). (The invariance with respect to Hq , i.e. to the geodesic flow, is proved in the same 
way as in the Riemannian case, and the invariance with respect to Vj is straightforward.) 

o 

Remark 3.2.3 (i) The vertical terms could be seen as an effect of the matter or the 
radiation present in the space-time M.. The S-diffusion ($ s ) reduces to the geodesic flow 
in the regions of the space where S vanishes, which happens in particular for empty space- 
times M in the cases S = -g 2 R(£), or S = g 2 £(£,£), or also S = -g 2 R(£) e K£ ^/ R ^ 
(for positive constant k) for example. 

(ii) As well as for the basic relativistic diffusion, the law of the S-relativistic diffusion 
is covariant with any isometry of (A4,g). The basic relativistic diffusion corresponds to 
S = g 2 > , and the geodesic flow to S = . 

(Hi) In [B] is considered a general model for relativistic diffusions, which may be 
covariant or not. Up to enlarge it slightly, by allowing the "rest frame" (denoted by z in 
[B]) to have space vectors of non-unit norm, this model includes the generic S-diffusion 
(compare the above equation (27) to (2. 5), (3. 3) in [B]). 

3.3 The ^-diffusion 

We assume here that the scalar curvature R = is everywhere non-positive on 

Ai, which is physically relevant (see [L-L]), and consider the particular case of Section 3.2 
corresponding to S = —g 2 R(^), with a constant positive parameter g. 

In this CclSG, clS its central term clearly vanishes, Equation (27) takes on the simple form : 

d 

d<5> s = H (<f> s ) ds + g J2V-R(®s)V j (<f> s ) o dwi . 

3=1 

3.4 The ^-diffusion 

We assume that the Weak Energy Condition (recall Section 2.2) holds (everywhere on 
T X A^), which is physically relevant (see [L-L], [H-E]), and consider the particular case of 
Section 3.2 corresponding to S = g 2 S = g 2 £(£,£) = g 2 T 00 . 

We call energy relativistic diffusion or ^-diffusion the G(.M) -valued diffusion process ($ s ) 
we get in this way, as well as its T 1 A1-valued projection 7Ti(<E» s ). 

The following is easily derived from Lemma 2.2.2 and Formula (10). As a consequence, the 
central drift term in Equation (27) is a function of the Ricci tensor alone when S is. 
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Lemma 3.4.1 We have VjR k = 5 0i R k - R\ + b\ Rij - 5 k R 0i , for < i, k < d and 
1 < j < d . In particular, VjR = , and V 3 £ = VjT 00 = V}i? o = 2R j . 

By Lemma 3.4.1, the drift term of Equation (30) which involves the derivatives of S 

d 

equals here: g 2 ^ Roj(£,s) e k (s) ds . As we have Rqj = R mn e™ e" by (16), we get the 
i=i 

d 

alternative expression : g 2 R 0j (Q ej(s) ds = g 2 R mn {Cs) f m € ~ ds ■ 

Another expression is got by using Einstein equation (17), or equivalently, by computing 

BE . . . 

directly from (30) and (12) : - g ki ] -r- = 2^ - g U \ f lm ^ = 2[E£- f where 

the notation (TC,) k = T^ m has the meaning of a matrix product. 
Hence, Formula (30) of Proposition 3.2.2 expressing d^ s reads here : 

dg = dM k - r*(&) £ H ds + ^S s 4 fc ds + g 2 R mn (is) t g C - 9 kn (Q] ds , (31) 
or equivalently : 

di = dM s - r. ■(£,) C g ^ + £? 2 (f + 1) ^ & 'Is - ,r' T & • (32) 

We can then compute the equation satisfied by the random energy S s . In particular, 
the drift term discussed above for d^ k contributes now for : 

2g 2 f km C[Si - fi} k = 2g 2 (£f km i m i k - f km C[f£] k ) = 2g 2 (S 2 - [T& [f£] k ). 
This leads to the following, to be compared with Corollary 3.1.2. 

Remark 3.4.2 The random energy S s := £(£ s ,£ s ) associated to the E- diffusion (<3> s ) 
satisfies the following equation (where V„ := v^Vj) : 

dE s = V^£(&, is) ds + (d + 2) g 2 E 2 ds - 2g 2 g(f£ s , ft) ds + 2g dM £ s , 

with the quadratic variation of its martingale part dM £ given by : 

[dS a , dS s ] = 4g 2 [dMf, dM £ s ] = 4g 2 [E 2 - g(f£ s , ft)] E s ds . 

Note that the diffusion coefficient [£ 2 — g(T£ s ,T£ s )] appearing in the quadratic variation 
[d£ s , d£ s ] is necessarily non- negative, which can of course be checked directly. 

Remark 3.4.3 Case of Einstein Lorentz manifolds. 

The Lorentz manifold M is said to be Einstein if its Ricci tensor is proportional to 
its metric tensor. Bianchi's contracted identities (see [H-E] or [K-N]), which entail the 
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conservation equations VkT jk = , force the proportionality coefficient p to be constant 
on M. . Hence : 

Rtm(£) — P 9£m{£,) , for any £ in M. and < £, m < d . 
Then the scalar curvature is R(£) — (d+ l)p , and by Einstein Equations (17) we have: 

TUO = (A - ^P)gem(0 =■ -pqUO ■ 

Hence Equation (19) holds, with q = : we are in a limiting case of perfect fluid. 
Moreover, R and 8 are constant, so that on an Einstein Lorentz manifold, the .R-diffusion 
and the ^-diffusion coincide with the basic relativistic diffusion (of Section 3.1). 

4 Warped (or skew) products 

Let us consider here a Lorentz manifold (Ai,g) having the warped product form: 
A4 — I x M , where I is an open interval of R + and (M, h) is a Riemannian manifold, is 
endowed with the Lorentzian pseudo-norm g given by : 

ds 2 := dt 2 - a(t) 2 \dx\ 2 h , (33) 

or equivalent ly 

9ok ■= S ok and := -a(t) 2 h i:j (x) , for < k < d, 1 < i, j < d. (34) 

Here £ = (t,x) E I x M denotes the generic point of M, and the expansion factor a is a 
positive C 2 function on I . The so-called Hubble function is : 

H{t) := a'{t)/a{t) . (35) 

This structure is considered in [B-E], which contains most of the following proposition. 

Proposition 4.1 Consider a Lorentz manifold (Ai,g) having the warped product form. 

(i) Rs curvature operator 1Z is given by : for u, v,w,a G C 2 (J) and X, Y, Z,A G T(TM), 

(n ((ud t + X) A (vd t + Y)) , (ad t + A)A (wd t + Z)} = -a 2 (/C (X A Y) , A A Z) 

+ aa" h(uY -vX,aZ - w A) + (aa) 2 [h(X, Z) h(Y, A) - h(X, A) h(Y, Z)\ , (36) 
/C denoting the curvature operator of (M, h) . 

(ii) Denoting by Ric the Ricci operator of (M, h) and by (■,■) the standard canonical inner 
product ofM. d , we have: 

a" 

(mcti(vdt + Y), wd t + Z) = (Ric(Y), Z) + [id - 1) \a'\ 2 + aa"} h(Y, Z)-d — vw. (37) 

' a 

(in) If V 9 and V h denote the Levi-Civita connections of(A4,g) and (M,h) respectively : 
V U+x) (.vdt + Y) = [uv' + a(t)a'(t) h(X, Y)} d t + W h x Y + H(t) (uY + v X) . (38) 
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Remark 4.2 In any warped local chart (t,x^) : for 1 < m,n,p,q < d (K mnpq denoting 
the curvature tensor of (M, h), and K mp its Ricci tensor) and < k < d , we have 

T 00 (g) = 0; T p 0n (g) = H 5? ; T° mn (g) = aa> h mn ; T? mn (g) = T? mn (h) , (39) 
K-onkg = 5 0k a(t)a"(t) h nq ; H mnpq = {aa') 2 (t) [h mq h np - h mp h nq ] - a 2 (t)K mnpq , (40) 



Rok = ~ S k d x a" (t) / a(t) ; R mp = K mp + [(d - 1) \a'\ +aa"]h 



mp ■ 



(41) 



Let us outline the proof for the convenience of the reader. 
Proof We get (38), and then (39), by using Koszul formula: 

2 h(V h x Y, Z) = Xh(Y, Z)+Yh(X, Z)-Zh{X, Y)+h([X, Y], Z)-h([X, Z], Y)+h([Z, Y],X), 
for both V h and V 9 . Hence, 

Ke t+ x)Ka t+ v)^ + Z ) = V U + x) (V + ««' h (Y, Z)]d t + V h Y Z + H [vZ + wY]) 



Therefore 



V {udt+xy V (vdt+Y) 



(wd t + Z) 



And since 



V! [(ud t+ XU vdt+ Y)]( wd t + Z )= ^[uv'-u'v]d t+ [X,Y])( Wdt + Z ) 

= ([uv' - u'v]w' + act h([X, Y],Z]^d t + V\ x;Y] Z + H {[uv' - u'v] Z + w[X,Y\), 
we get 



V («9t+X)' V (vdt+Y) 



(ud t +X),(vdt+Y)] 



(wd t + Z) 



= aa"h(uY-vX, Z)d t +([V h x , V Y }-V h [XtY] )Z+—w{uY-vX)+a l2 [h{Y, Z)X-h(X, Z)Y\. 
By Formula (9), this entails Formula (36), which is equivalent to (40), by (15). 
Then, denoting by (e± : . . . , e^) an orthonormal basis of (M, h) : 

which yields (37) ; which is in turn equivalent to (41), by (16). o 

Corollary 4.3 A Lorentz manifold (Ai,g) having the warped product form is of perfect 
fluid type (recall Section 2.4) if and only if the Ricci operator of its Riemannian factor M 
is conformal to the identical map : Ric = Q x Id , for some Q G C°(M). 

If this holds, we must have : U = d t , p + q = —d (a" /a), and q = fl a~ 2 — (d — 1)H' . 
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Proof By (21) and (37), this happens if and only if for any v E C°(I),Y E T(TM) : 

(Ric(Y), Y) + [(d - l)a' 2 + aa"}h(Y, Y) - d(a"/a)v 2 = q g(U, vd t + Y) 2 +p [v 2 - a 2 h(Y, Y)] , 

which (seen as a polynomial in v) forces U — d t , and then splits into p + q = —d (a" /a) 
and (Ric(Y), Y) = —[(d — l)a' 2 + aa" + pa 2 ] h(Y,Y). This latter equation is equivalent 
to p — —[(d— 1)H 2 + (a" /a) + ila~ 2 ] , and then, using p + q = —d(a"/a), to q = 
tta~ 2 - (d- 1)H'. o 



Corollary 4.4 Consider a Lorentz manifold (A4,g) having the warped product form. The 
energy (12) at £ = (i,x) E T^M = T^ x) (I x M) equals : 



£(t, £) = (Ric(x),x) - (d - 1) H'(t) (i 2 -l) + \d{d- 1) H 2 (t) + 



R 



M 



2a 2 (t) 



(42) 



Then the weak energy condition is equivalent to the following lower bounds for the Ricci 
operator and the scalar curvature R M of the Riemannian factor (M, h) : 



inf > ( d _ i) sup { a a " - a ' 2 }; R M > - d (d - 1) inf a' 2 . 

wgtm h(w,w) ~~ ' i ~~ i 

And the scalar curvature R of (Ai,g) equals: 



R = - a - 2 R M - d 



(d-1) |^| 2 + 2^ 

V /la 1 a 



(43) 



(44) 



so that its non-positivity is equivalent to the lower bound on the scalar curvature R M : 

R M > -dx mi{(d-l)a' 2 + 2aa"} . (45) 



Proof From (37), we compute the scalar curvature of M. : 



R = (Ricci(^), dt) v - a~ 2 ^(Ricci(e 3 -), ^\ = -a~ 2 R M ~ dUd - 1) |£ \ 2 + 2 

3=1 

On the other hand, by (37) we have at £ = (t, x) E T£M = T^ x) (I x M) : 
(Ricci(0,0, = (Ric(x),x) + [(d - 1) \af + aa"} ^ - d^i 2 

= (Ric(x) , x) - (d - 1) H'(t) (i 2 - 1) - d (t) . 

Thence Formula (42). Then, the weak energy condition holds if and only if for any t E I 
and w E TM : 



(Ric(w),w) > (d - 1) a 2 {t)H'{t) h(w, w)-\d{d-l) H 2 {t) - 



M 



R 

2aHt) 
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By homogeneity with respect to w , this can be split into the following lower bound for 
the Ricci operator of the Riemannian factor M : 

„ (Ric(w), w) , , . , „ ,0-. 
inf v ,/ J \ 1 > (d- 1 ) sup {a a" -a' 2 }, 

w&TM h(w,W) I 

together with the condition particularised to w — , which yields the following lower 
bound for the scalar curvature of the Riemannian factor M : 

d(d - 1) H 2 (t) + a(t)~ 2 R M > , or R M > - d (d - 1) inf a' 2 . o 

5 Example of Robertson- Walker (R-W) manifolds 

These important manifolds are particular cases of warped product : they can be written 
M = /xM, where / is an open interval of M + and M G {§ 3 , M 3 , H 3 }, with spherical co- 
ordinates £ = (t, r, ip, ip) (which are global in the case of M 3 ,H 3 , and are defined separately 
on two hemispheres in the case of S 3 ), and are endowed with the pseudo-norm: 

9(1 ■= i 2 - a{tf (y^t 2 + r V + r 2 sin 2 ^ ^ , (46) 

where the constant scalar spatial curvature k belongs to {—1,0,1} (note that r G [0,1] 
for k = 1 and r G R+ for k — 0, — 1), and the expansion factor a is as in the previous 
section 4. Note that we have necessarily t > 1 everywhere on T l Ai . 

By (36), we have the curvature operator given by: 

K {{ud t + X) A (vdt + Y)) , (ad t + A) A (wd t + Z)\ 

= aa" h(uY - vX, aZ - wA) - a 2 (a' 2 + k) [h(X, A)h(Y, Z) - h(X, Z)h{Y, A)} . 
By (37), the Ricci tensor ((Rem)) is diagonal, with diagonal entries: 

a"J^_ A(t) 
t(t) 

and the scalar curvature is R — —6 [cx(t) a"(t) + a'(t) 2 + k] a(t)~ 2 . 

The Einstein energy-momentum tensor Re m — \Rge m = Te m is diagonal as well, with 
diagonal entries : 

/ a'{tf + k ^ -A(Q ) _ A ^ r 2 ; r 2 gin 2 \ ) with A ^ ._ 2a (t)a"(t) + a\tf + k . 

Hence, we have 

ftm - a{ty 2 A{t) ge m = 2 [ka(ty 2 - H'(t)\ l { e= m = } ■ 

18 



(- 3 , - A{t ) - , A(t) r 2 , A(t) r 2 sm 2 V ) , where A(t) := a(t) a"(t) + 2 a'{t) 2 + 2A; , 
V a(t) 1 — kr 2 J 



Thus, in accordance with Corollary 4.3, we have here an example of perfect fluid : Equation 
(19) holds, with 

Uj = q, -p(0 = ka(ty 2 + 2H'(t)+3H 2 (t) , q(g) = 2 {ka(t)~ 2 - H'(t)] , (47) 
p(0 = -2 [2k a{t)- 2 + H\t) + 3H 2 (t)]/(d - 1) . 

Note that 

As = CiUC = i s and E s = 2 [ka(t s )- 2 - H'(t s )} i 2 s - p(Q . (48) 

By Corollary 4.4 (or by Remark 2.4. l(ii) as well), the weak energy condition is equivalent 
to : a' 2 + k > (aa") + . 

We shall consider only eternal Robertson- Walker space-times, which have their future- 
directed half-geodesics complete. This amounts to I = M* + , together with 



a 



in such Robertson- Walker model) , we have in particular : 



oo . In the case of the basic relativistic diffusion (solving Equation (23) 



dt s = Q\J t 2 - 1 dw s + 1 8 ds - H{t s ) [t 2 s -l\ds. (49) 

5.1 S-relativistic diffusions in an Einstein-de Sitter-like manifold 

We consider henceforth the particular case I — ]0, oo[ , k — , and a(t) = t c , 

with exponent c > . Note that such expansion functions a can be obtained by solving 
a proportionality relation between p and q (see [H-E] or [L-L]). 

Thus q = 2ct~ 2 , p = (2-3c)ct- 2 , R = -6 c (2c - 1) r 2 , E = ct~ 2 (2 i 2 + 3c - 2) . 

Note that the weak energy condition holds. The scalar curvature is non-positive if and 
only if c > 1/2 , and the pressure p is non-negative if and only if c < 2/3 . 

Note that the particular case c = | corresponds to a vanishing pressure p , and is 
precisely known as that of Einstein-de Sitter universe (see for example [H-E]). And the 
analysis of [L-L] shows up precisely both limiting cases c = | and c—\. 

5.1.1 Basic relativistic diffusion in an Einstein-de Sitter-like manifold 

In order to compare with the other relativistic diffusions, we mention first for the basic 
relativistic diffusion (of Section 3.1), the stochastic differential equations satisfied by the 
main coordinates i s and r s , appearing in the 4-dimensional sub-diffusion (t s , i s , r s , r s ). By 
(49), we have, for independent standard real Brownian motions w,w: 

di s = Qsji 2 s - 1 dw s + & i 8 ds - — (i 2 s -l)ds; (50) 
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Almost surely (see [A]), lim i s — oo , and x s /r s ~ x s /\x a \ converges in § 2 . 

5.1.2 /^-diffusion in an Einstein-de Sitter-like manifold 

With the above, Section 3.3 reads here, for the i?-relativistic diffusion, when c > 1/2 : 



di s = Q dM s + 9c (2c - l)g 2 t~ s 2 £s ds - V^Q £ £ rfs , 



(52) 



with the quadratic covariation matrix of the martingale part dM s given by : 

Q- 2 [dg, dg] = 6c (2c - 1) [g g - /%)] t; 2 (fa , for < k, £ < d 
In particular, we have for independent standard real Brownian motions w, w : 

9^ 2 c(2c- 1) . 



dt s = -\/Qc{2c- l){i 2 s - l)dw s 4 



tl 



t s ds (t s — 1) ds ; 

t.q 



df x = 



+ 



Qy/6c(2c-l) 

9^ 2 c(2c- 1) 



dtVc 



(53) 
(54) 



tl 



r s ds + 



i 2 - 1 

t 2c 



ds 2c • 
t* r, ds . 



As the scalar curvature R s = 6c (1 — 2c)/t 2 vanishes asymptotically, we expect that almost 
surely the /^-diffusion behaves eventually as a timelike geodesic, and in particular that 
lim t* — 1 . 



5.1.3 ^-diffusion in an Einstein-de Sitter-like manifold 

Similarly, using (32), (47), (48), we have here £ f - T£ = 2(0 - H')(i 2 i -W) , so that 
Section 3.4 reads here, for the ^-diffusion : 

di a = gdM s + 3 -^ tf (2 i 2 s + 3c - 2) g <fa + 2f? 2 c t; 2 (i s g - C/ s ) t a ds - (&) g g <fa , (55) 

with the quadratic covariation matrix of the martingale part dM s given by : 

Q- 2 [dil dg] = c [g g - «^(£ s )] (2 i 2 + 3c - 2) t; 2 da , for < k,£ < d . 

In particular, we have for some standard real Brownian motion w : 



di, 



(2t2-2 + 3c)(t2-l)d Ws + c 
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V(t 2 -i + !)S- 

s 



10 y « t. 



(is ; (56) 
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dr s = ^\l2ti -2 3r 



1 



dw s + a / — 

V 7 ^! V t 2c tl - I 



+ £ 2 c (5 i 2 s - 3 + f ) ^ - ^ t s r s 4 



2c , 



^-1 

* 2c 



dw s 
ds 

Ts 



(57) 



Remark 5.1.4 Comparison of ^-diffusions in an E.-d.S.-like manifold 

Along the preceding sections 5.1.1, 5.1.2, 5.1.3, we specified to an Einstein-de Sitter-like 
manifold the various H-diffusions we considered successively in Sections 3.1, 3.3, 3.4. Re- 
stricting to the only equation relating to the hyperbolic angle A s = t s , or in other words, 
to the simplest sub-diffusion (t s ,t s ), this yields Equations (50), (53), (56) respectively. We 
observe that even in this simple case, all these covariant relativistic diffusions differ notably, 
having pairwise distinct minimal sub-diffusions (with 3 non-proportional diffusion factors). 



5.2 Asymptotic behavior of the /^-diffusion in an 
Einstein-de Sitter-like manifold 

We present here the asymptotic study of the i?-diffusion of an Einstein-de Sitter-like 
manifold (recall Sections 5.1, 5.1.2). We will focus our attention on the simplest sub- 
diffusion (t s ,i s ), and on the space component x s G IR 3 . Recall from (48) that i s = A s 
equals the hyperbolic angle, measuring the gap between the ambient fluid and the velocity 
of the diffusing particle. Recall also that, by the unit pseudo-norm relation, i s controls 
the behavior of the whole velocity £ s . We get as a consequence the asymptotic behavior 
of the energy £ s . As quoted in Section 5.1.2, we must have here c>\- 

Note that for c = \ , the scalar curvature vanishes, and the i?-diffusion reduces to the 
geodesic flow, whose equations are easily solved and whose time coordinate satisfies (for 
constants a and s ) : 

s — s = a/^s (t s + a 2 ) — a 2 \og[y/t^ + \/t s + a 2 ] , whence t s ~ s . 

The proofs of this section (and of the following one) will use several times the elementary 
fact that almost surely a continuous local martingale cannot go to infinity. 

The following confirms a conjecture stated at the end of Section 5.1.2. 
Proposition 5.2.1 The process i s goes almost surely to 1, and E s — > , as s — > oo . 
Proof By Equation (53), 

log ^ - 3, 2 c(2c - 1) f (2 + r 2 )% + c f (1 - if) l f dr 
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H ■ dr 

is a continuous martingale with quadratic variation 6g 2 c(2c — 1) / (1 — i~ 2 ) — ■ Hence, 

Ji t T 

since t T > 1 and therefore t T > r , the non-negative process 

logt s + c / {l-i-^dr 
Ji t T 

converges almost surely as s — > oo . This forces the almost sure convergence of the 
f°° i T 

integral : / (1 — i~ 2 ) — dr < oo , and of i s , towards some G [1, oo[ . This implies in 

Jl tr 

f°° ■ dr 

turn t T = 0(t), hence / (1 — i~ 2 ) — < oo , whence finally = 1 . o 

Ji r 

Consider now the functional a := t c , which is constant along any geodesic. 

Lemma 5.2.2 For c > \, the process a s := t c s \Ji 2 s — 1 goes almost surely to infinity, and 



3C 



ds 

cannot vanish. Moreover, for any e > we have almost surely : I tf~ 2 —— — < oo 



at 



Proof We get from Equation (53) : 

da s = - ^Qc (2c - 1) (aj + tf) dw s + 3g 2 c (2c - 1) 3 °% + 2 ^ ds , 
t s t s a s 

and then for any e G ]0, 1] and for some continuous local martingale M : 

< ar = ar £ ~M S - 3eg 2 c (2c - 1) f [2 ~ " ] a \ + }\ - " ] ^ dr. 

The signs in this last formula, and the fact that almost surely a continuous local martingale 
cannot go to infinity, imply the convergence of the last integral and of the martingale term 
M s , entailing the almost sure existence of a finite limit , hence of G ]0, oo], and the 

d: 



almost sure convergence of the integral / — ^- — ^— ^- < oo . Now, by Proposition 5.2.1, 

ai +£ t±- 2c 



f°° dr [°° dr 
this implies / — - — < / — — — — — < oo , hence = oo . Finally, the equation for 

a~ e forbids also the existence of a finite zero s for a s . Indeed, s /* s would force the 
martingale term of this equation to go to — oo, which is impossible, o 

The following reveals the asymptotic behavior of the space component (x s ) for c > | . 
Proposition 5.2.3 For c > \, the space component converges almost surely (as s — > oo) : 
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Proof (i) Let us consider the non-negative process u s := t s (i 2 — 1), which is constant 
for c = \ , and cannot vanish for c > \ , by Lemma 5.2.2. By Equation (53), 



u. 



+ (2c-l) f S ^dr- 6g 2 c(2c-l) f [4t 2 -l] — 



is a continuous local martingale. Then, for any e > , 

s -6g 2 c[2c-l] 1 \ > + [2c-l + e 1 y 



" POO 













+e J / +l+e L J / +l+e 

Jl L T Jl L T 

is a continuous local martingale. By Proposition 5.2.1, the central term converges almost 
surely. This implies that J < oo and that — converges, almost surely. 

(ii) By the unit pseudo-norm relation, we have t 2 s c \x s \ 2 = t 2 — 1 = u s /t s . Let us apply 
(i) above with e :— c — \ , to get : 

"CO POO POO 

This proves that x s = x\ + / x T rfr — )■ X\ + / i s ds£l , almost surely as s — >■ oo . o 
In the case c = \ of the /^-diffusion being the geodesic flow, we have 

r s = a/6 2 /a 2 + (a + o(l)) logs ~ a/o log s as s — >■ oo , 

which shows that Proposition 5.2.3 does not hold for the limiting case c = \ ■ 

To compare the /^-diffusion with geodesies, note that (as is easily seen ; see for ex- 

ample fAl) along any timelike geodesic, we have x s = X\ + - — r / — — (and = 

' Fi| Ji tf v l Xs l 

lHj) , which converges precisely for c > \ ; and along any lightlike geodesic, we have 
x s = rci + / — ~ V x si+^ (and p^, = tPt), which converges only for c > 1 . 

On the other hand, for c < 1 , the behavior of the basic relativistic diffusion proves to 
satisfy (see [A]) : r s ~ / — " 2e — ► oo (exponentially fast, at least for c < 1). 

Hence, the /^-diffusion behaves asymptotically more like a (timelike) geodesic than like 
the basic relativistic diffusion. However, owing to Lemma 5.2.2, the asymptotic behavior 
of the /^-diffusion seems to be somehow intermediate between those of the geodesic flow 
and of the basic relativistic diffusion. 
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5.3 Asymptotic energy of the ^-diffusion in an E.-d.S. manifold 



We consider here the case of Section 5.1.3, dealing with the energy diffusion in an 
Einstein-de Sitter-like manifold, and more precisely, with its absolute-time minimal sub- 
diffusion (t s ,t s ) satisfying Equation (56), and with the resulting random energy: 

S s = ct- 2 (2 t 2 s + 3c - 2) = 2c {ts/tsf + 0(s- 2 ). 
Let us denote by ( the explosion time : ( := sup{s > | t s < oo} G ]0, oo] . 

Lemma 5.3.1 We have almost surely : either lim t s — 1 and ( = oo , or lim i s — oo . 

s— s-*C 

Proof By Equation (56), 

1 f sA( > dr /" sAC dt 

~r~ c J [^'Ut + sc [ 3 + W + 2 W 

tsA( J S l T J SO l T 

f sAC _ dr 

is a continuous martingale with quadratic variation 2^ 2 c / [l + [l — p-] — ■ Hence, 

since i T > 1 and therefore t T > r , the process 

"SAC 



Js l T 



converges almost surely as s — > oo . As < t sA( * < 1 , this forces the almost sure 

f> dr 

convergence of the integral : / [l — — < oo , and of t S AC > towards some G [1, oo]. 



Moreover, the convergence of the integral forces either ( < oo and then = oo , or 
C = oo and then G {1, oo}. o 

The asymptotic behavior can, with positive probability, be partly opposite to that of 
the preceding /^-diffusion : 

Proposition 5.3.2 From any starting point (t S0 ,t SQ ), there is a positive probability that 
both A s = t s and the energy S s explode. This happens with arbitrary large probability, 
starting with t So /t SQ sufficiently large and to bounded away from zero. 

On the other hand, there is also a positive probability that the hyperbolic angle A s = t s does 
not explode and goes to 1, and then that the random energy S s goes to 0. This happens 
actually with arbitrary large probability, starting with sufficiently large t S0 /i SQ . 

Proof Let us set A s := t s /i s > . From the above proof of Lemma 5.3.1, we get directly : 
X SA( - X so = M sAC + [1 + c] / [l-^ m ]dr- Q 2 c [3 + 1^ + ^]— (58) 

J So T J SO T T 
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/•SAC 

< M sK + [1 + c](s A C - s ) - g 2 c / 

J SO 



- 3 l + 2i~ 2 ] 

— , dr, (59) 

(M s ) denoting a martingale having quadratic variation 2g 2 c / [l + ^f^][l — i?\dr. 

J So 

(i) Let us first start the time sub-diffusion from (t So ,i SQ ) such that t SQ > So > 1 and 
t so >nm t SQ , with fixed n > 2 and m > 2 + , and consider 

T := £ A inf {s > s | m \ s > 1}. Thus, we have on [s , T] : A" 1 > m , and then 

3 — 1+2 t T > g m _ 3 > 3 _|_ i+c _ Therefore, by Inequality (59), we have almost surely for 
any s > s : 

< A sAT < A so - 3g 2 c (s A T - s ) + M sAT . 
Integrating this inequality and letting s / oo yields : 

i ¥ (t S0 ,ijT <(]< liminf E[A SAT ] < A So < ^ , whence F^JT = (} > 1 - \ . 
Moreover, almost surely on the event {T = (}, the above inequality 

< Kac < K ~ 3f? 2 c (s A C - so) + M sK 

implies clearly (using that a continuous martingale almost surely cannot go to infinity) 
( < oo , and by the previous lemma that = oo . Then (58) implies the convergence of 
A sA £ to some A^ G K + . 

Furthermore, = oo and A^ > for finite ( would imply trivially — oo , whence its 
logarithmic derivative should explode, which leads to a contradiction. 

This proves that we have P(t so ,t so ) [C < 00 > \ = 0] — 1 — 77,-1 ■ 

Since (by the support theorem of Stroock and Varadhan, see for example Theorem 8.1 
in [I-W]) from any starting point the sub-diffusion (t s ,t s ) hits with a positive probability 
some (t so ,t so ) as above, we find there is always a positive probability that i s and £ s 
explode (together). 

(ii) Let us now start the time sub-diffusion from (t Sl ,t Sl ) such that nm't Sl < t Sl , 
with fixed m' > 9g 2 c (1 + c) and n > 2 + c , and consider T' := ( A inf {s > s± | \ s < m'}. 

By Equation (56) we have at once : almost surely, for any s > s 1 , 

psAT' psAT' ns AT' 

Kat> = K! + 5g 2 c / [1 + ] A; 3 dr-[l + c] A; 2 dr + c r 2 dr + M' AT , , 

J Sl J S\ J Si 

(M' s ) denoting a martingale having quadratic variation 2g 2 c / [l + 3 f^][l — dr. 

J Sl 

Since on [s u V] we have A; 1 < 1/rri, and then 5^ 2 c[l + ^f] X; 1 < 5g2c ^ +c) < 5/9 , we 
get: 

< A S - A V < A," 1 - c / A; 2 dr + — + M' SAT , . 

J si r si 
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f T ' 

/ A r 2 dr < oo and X S ^ T , — > X T } G 



(which implies moreover \ T } = 



This entails 

almost surely on {T' = 00} ), and 

h < CI < E (t , lAl) [A^ 1 ] < A" 1 + ^ < ^ • 

Hence, we get P^ t t = C] > 1 — ^ > • Furthermore, as in (i) above, tr 1 = 00 
and Xt' > for finite T' is impossible, which excludes T' = ( < 00 . Therefore 

P(t n ,t n) r = C = oo] >l-3±s >0. 

Then from the equation for A , using that [3 + f ]g 2 c\- 1 < {&+ ^f c < \ on [si,T'], we 
get almost surely : 

r sAT " dr 

A s at' - A S1 > (s A T' - Sl ) -[3 + ^]g 2 c — + M sAT , > |(s A T' — si) + M sAT , , 

which shows (since [M s , M s ] = 0(s)) that almost surely {T' = ( = 00} C {\ s — > 00}. On 
this same event, by Equation (56) we have almost surely for any s > s±: 



t a - <si = Q 



2c[l + ^][l-±]^dw T + c 



A T 



5g 



1 1 9c-10 
2 10*2 



l-t 



1-2 



A T 



which shows that i s cannot go to infinity, since this would forbid the last integral, and 
then the right hand side, to go to +00 . Hence, by Lemma 5.3.1, we obtain that almost 
surely {T' = ( = 00} C {i s — > 1}. The proof is ended as in (i) above, by applying the 
support theorem of Stroock and Varadhan, and by taking n arbitrary large, o 



6 Sectional relativistic diffusion 

We turn now our attention towards a different class of intrinsic relativistic generators 
on G(A4), whose expressions derive directly from the commutation relations of Section 2.2, 
on canonical vector fields of TG(A4). They all project on the unit tangent bundle T X M. 
onto a unique relativistic generator T-i l curv , whose expression involves the curvature tensor. 
Semi-ellipticity of T-L\ urv requires the assumption of non- negativity of timelike sectional 
curvatures. Note that in general T-L l curv does not induce the geodesic flow in an empty 
space. 

6.1 Intrinsic relativistic generators on G(A4) 

We shall actually consider among these generators, those which are invariant under the 
action of SO(d) on G(Ai). To this aim, we introduce the following dual vertical vector 
fields, by lifting indexes : V ij := rf m rf n V mn . Note that V j = V 0j = -V 0j = -Vj , and 
that V* j = for 1 < i,j < d. We consider again a positive parameter g. 
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Proposition 6.1.1 The following four S 0(d) -invariant differential operators define the 
same operator 1-L l curv on T l M : 

d d 

# - 4 E ([#0' H A yj + H A) ; H » + f E^ ; 

j=l l<j,k<d l<i,3<d 

Note that {T-i l curv — C°) is self-adjoint with respect to the Liouville measure ofT x M.. 

The proof will be broken in several lemmas. We begin with the following general and useful 
computation rules, derived from Sections 2.1 and 2.2. 

Lemma 6.1.2 For < j, k, £ < d , we have : 

V l TZ l3 M = 5 k R] - 5* R) + (1 - d)K 0j u + K e j0 k - 1Z k j0 e ; 

V*n 0l ke = 5 £ R k -5 k R e ; [[H^H&V*] = (d - 1)[H , H 3 ] . 

Proof We get the first formula by multiplying by i] ip the formula of Lemma 2.2.2, and 
particularising to q = . As to the second one, by particularising the latter to j = and 
changing sign, we get : 

■friTj ki r£ ryk rfc ryi i Tjk I -jjl k 
V 'MM — OQ-tiQ — K + /t oo — K> 00 • 

Then, we note that 7^. £ o fc = 7^0% = 7^ fc oo £ - Finally, the last formula derives from the 
second one and from the commutation relations (4) and (6), as follows: 

[ [Hi, h 3 ], v l ] = \ \n l3 M v ke , v l ] = \ n i3 u [v ki , v l ] - \ (v l n l3 ke )v ke 

= \ K i3 u (Sjy t - 5}V k + V ip ( m V pk - r^)) - (6 k R] + ^ TZ, 3 M + 7^ fe ) V u 

= ^.Ve+n^Vpk-E^VtH^^Vkt-n^VM = (^)n 0j ke v ke = (d-i)[H ,H 3 }. o 

We get then first the following. 
Lemma 6.1.3 On C 2 (T 1 M), we have [[H ,H 3 },V j ] = 0, and 

d d d 

3=1 J 3=1 3=1 l<j,k<d 
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Proof Using the commutation relations (4) and (6), that Vu = on C 2 (T l M) for 
1 < k, £ < d , and (10), we have on one hand : 

[H ,H 3 }V = \n 0j u v H v j = lK 0j H ([V ke ,V j ] + V j V ke ) 
= \ K 0j H V ij (VikV e - ViiV k + r) 0e V ik - Vok V t£ ) + n 0j oe V>V e 

= -n 0j kj v k + n ih v ik + n 0j ok v j v k = n 0j ok v j v k -R k v k . 

On the other hand, using this first part of proof and Lemma 6.1.2, we get : 

[[H^Hjiv*] = \[n Gj u v u ,vi] = ±n 0j M [v ke ,v>]-±(vm 0j ke )v ke 

= -R k V k - \{8i R k - 5 k R e )V ke = -R k V k + R k V k =0. 
Using [H , Hj]V j + V j [H , H 3 ] = 2 [H , Hj]V j - [ [H , H 3 ], V j ] ends the proof, o 

We get then the following. 
Lemma 6.1.4 On C 2 (T l M), we have £ [[H h H 3 ], V ij ] = , and 

l<i,j<d 

£ {[Hi, Hj] V« + V* [H U Hj]) = - 4 n° curv . 

l<i,j<d 

Proof As for the proof of Lemma 6.1.3, we use the commutation relations (4) and (6), 
that V u = on C 2 {T l M) for 1 < k,£ < d, (10), Lemmas 2.2.2 and 6.1.2, and the 
symmetries of the Riemann tensor. We have thus on one hand and on C 2 (T 1 A^) : 

[[H p ,H 3 },v ij ] = \[n pj kl v M ,v mn ]rf m rf n = \n P3 M [v M ,v mn ]^ n - \{v mn ii V3 kl )if m ^ n Vu 



= \ (TZ p m£ V en + R k p V k l - K p nki V kn +K e p V e i ) + l[(d+ l)1l\ ke + if k R l p + TZ p m - V u R k - TZ p ak ] V kl 
= §[- Tl kU + r] ie R k - TZ p eki - if k R l p + (d + l)K iki + r] tk R e p + K kii - if 1 R k - TZ p eik \ V M 

In particular, we get [ [H h Hj], V*\ = (^) K\ M V u = . 
And on the other hand, on C 2 (T 1 .M) again: 

[H t ,H,]v^ = Iff^yMK™, = rf^ n n l3 M v M v n = KJ M V M V 3 

= n jk \[V ke , V 3 ] + Vj V u ) = K jk \vjkVe - Vj eV k + m V jk - Vok V je ) + TZ jke V 3 V u 
= -2 R k V k + 2 K jk V jk + 2 TZ jM VjVt = 2 {TZ^ k VjV k - R k V k ) = -2 H° curv . o 

The final assertion relating to the Liouville measure is proved as in Theorem 3.2.1. 
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Proposition 6.1.5 In local coordinates, the second order operator T-L\ urv defined on T X M. 
by Proposition 6.1.1 reads: 

Hlurv = ee 4 4- + 4 cR k n 4r - i n p k 1)1 



^ SS ^3 Q^ k 1 2S -mnpg^y y ^ fc s y y 

Proof By Section 2.3, we have on C^^M) : 

8 2 8 8 2 8 8 2 8 

tttt n £ c n n I n n i z n 

VjVk = e ^ k deWo + ' ke °^ + e ° e °M^i +ej de V k = ejefc 5^ + ^ e °^' 
whence by Lemma 6.1.3 : 

= + K e > lk (including now 3 = 0) • 

On the other hand, by Formula (18) we have : 

= n 0abc V a 'v b °V ck = C e n a n mnp £ el e[ rj aj rf V ck = e™ e p e r a n mrp£ e[ rj aj rf\ 

whence 

Cj e q k = Co e$ e a lZ mr p£ e c r/ J r\ e q k = e e^TZmrpig g q = e eQ7?. m 
And in a similar way, by (16) : 

pi „« _ p m q p ij n _ m fm _ m -fi lq pn 

This and (13), (14) yield the wanted formula, whose coefficients depend only on (£,£) G 
T 1 ^, as it must be by the 50(d)-invariance underlined in Proposition 6.1.1. o 



6.2 Sign condition on timelike sectional curvatures 

The generator %\ urv defined on T l A4 by Proposition 6.1.1 is covariant with any Lorentz 
isometry of (Ai,g). Hence, it is a candidate to generate a covariant "sectional" relativistic 
diffusion on T 1 ^, provided it be semi-elliptic. 

As a consequence of Section 6.1, the intrinsic sectional generator we are led to consider 
on T X M is U\ urv , restriction of H + V j - R j0k VjV k ) . 
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Now, a necessary and sufficient condition, in order that such an operator be the gener- 
ator of a well-defined diffusion, is that it be subelliptic. 

We are thus led to consider the following negativity condition on the curvature : 

(lZ(u A v), u A v) < , for any timelike u and any spacelike v . (60) 

This condition is equivalent to the lower bound on sectional curvatures of timelike planes 

Ru + Rv: 

(TZ(u Av),u Av} r , > 
g(uAv,uAv) ~ 

since g(u A v, u A v ) :— g(u, u)g(v, v) — g(u, v) 2 < for such planes. 

Note that Sectional Curvature has proved to be a natural tool in Lorentzian geometry, see 
for example [H], [H-R]. 

We test this negativity condition on warped products, in Corollary 6.2.1 below. When 
this negativity condition is fulfilled, we call the resulting covariant diffusion on T X A^, which 
has generator Hl urv given by Propositions 6.1.1, 6.1.5, the sectional relativistic diffusion . 



Corollary 6.2.1 Consider a Lorentz manifold (Ai,g) having the warped product form. 
Then the sign condition (60) is equivalent to: a" < on I, together with the following 
lower bound on sectional curvatures of the Riemannian factor (M, h) : 

UC(X AY) ,X AY) . „ , 2l 

*,^« MX, X) h{Y\ Y) - h(X,Y)'' £ T ( " " } ' (W) 

Proof Let us denote by SR(U AV) the sectional curvature of the timelike plane associated 
with U AV. Recall from Section 6.2 that the negativity condition (60) reads simply 
SR(U A V) > , for any timelike U and spacelike V. By choosing a pseudo-orthonormal 
basis of such given timelike plane, we can moreover restrict to g(U, U) = 1 = — g(V, V) and 
g(U,V) =0. Setting U = ud t + X and V = vd t + Y , with u,v G C°(I) and X, Y e TM, 
we can thus suppose that : 

u = y/a 2 h(X,X) + 1 ; v = ^ a 2 h(Y, Y) - 1 ; uv = a 2 h(X, Y) , 

which implies a 2 h(X AY,X AY) + h(Y, Y) — h(X, X) = a" 2 , and 

h(uX - vY, uX - vY) = a 2 h(X AY,X AY) + oT 2 . 

Recall that h(X A Y, X A Y) := h(X, X)h(Y, Y) - h(X, Y) 2 . Now, by (36), this entails : 

SR(UAV) = a 2 (JC (X A Y) , X A Y)-aa"[a 2 h(XAY, X AY)+a~ 2 ] + \aa'\ 2 h(XAY, XAY) 

= a 2 \{IC(X AY),X AY}- (a a" — a' 2 ) h(X AY,XA F)j - a" /a. 
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Reciprocally, given any X, Y G TM such that 

a 2 h(X A F, X A F) + /i(F, F) - h(X, X) = a' 2 and h(Y, Y) > a~ 2 , setting 

u := v/a 2 /i(X,X) + 1 and v := ±v/a 2 /i(F,F) - 1, we get {uvf = a 4 h(X,Y) 2 , and 
thence a basis (U, V) as above. 

Hence, the negativity condition (60) is equivalent to : 

< a 2 ( (/C (X A Y) , X A F) - (a a" - a' 2 ) h(X A Y, X A F)) - a"/a , 

for any X,Y e TM such that a 2 /i(X A F,X A F) + /i(F,F) - /i(X,X) = a" 2 and 
/i(F, F) > a' 2 . 

Distinguishing between collinear and non-collinear pairs X, F, and denoting in the 
latter case by SK(X AY) the sectional curvature of the plane associated with X AY, this 
condition splits into both : (a"/ a) < , together with : 



a 2 h(X A F, X A F) 



< SK(X A Y) - (aa" -a' 2 ), 



for any non-collinear X, F G TM such that a 2 /i(X A F, X A F) + h(Y, Y) - h(X, X) = a~ 2 
and h(Y, Y) > a' 2 . We shall have proved that this condition is equivalent to the wanted 
inequality (61), if we show now that for any given a > , any given plane P in TM 
is generated by pairs X, F such that a 2 h(X A F, X A F) + h(Y, Y) - h(X, X) = a~ 2 , 
h(Y, Y) > a~ 2 , and such that h(X A F, X A F) is arbitrary large. 

Now, starting from an arbitrary {X , F } generating P, take F := — . y ° and 

ay/h(Y ,Yo) 



X:=q 



Y h(X ,Yo) V 



Then a 2 /i(XAF,XAF) + /i(F,F)-/i(X,X) = a~ 2 = h(Y, F), 



and h(X A F, X A F) = q 2 M ™g ro> is indeed arbitrary large, for arbitrary large q . 



o 



In particular, in an Einstein-de Sitter-like manifold, the sign condition (60) holds if and 
only if a" < , i.e. if and only if c < 1 . 

6.3 Sectional diffusion in an Einstein-de Sitter- like manifold 

We must have here c < 1 . By (40), we have : for < k < d and 1 < m,n,p, q < d, 



T^-Onkq ^Ofc C (c 1 ) t h n q and C t \hmq h n p h m p h n( ^\ t K 



mnpq • 



Using Cartesian coordinates (x^) instead of spherical coordinates (r, tp, vp) for the Euclidian 
factor M = M 3 , we have merely : 

T^-Onkq C (c 1) t (5ofe 5 n q and TZ-mnpq C t [5mq $np $mp $nq\ ■ 

Thence, for < k < d and 1 < m,n,p,q < d: 

K n k « = c(c-l) r 2c ~ 2 5 0k 5 n « and K m n p « = c(c-l) t 2c ~ 2 5 mp 5 tn 5 l<1 +c 2 1~ 2 [6« m 5 r ;-5 mp 5 n % 
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And Rij = Sij (3c (1 - c) t~ 2 5 0j + c (3c - 1) t 2c ~ 2 [l - 6 0j \), whence 

R k n = -5 k n cr 2 (3(c - l)8 0n + (3c - 1)[1 - Son}) ■ 
And by (39), the non- vanishing Christoffel coefficients are: 

r* = r% = c r 1 5* , and r°. = c t 20 ' 1 5 lj , for 1 < i, j, k < d , 
Therefore, by Proposition 6.1.5, we have: 



d_ 



9 



<9 2 



_ ij_d c . 2 _ d_ _ 2c • ■ <9 

" 4 t [ >dt T dxl 



2t 2[C L)t di 2t 2[6 ° L)X dxi 



g 2 c{c-l) 



2t 2 



(t 2 - 1) 



di 2 



g 2 c(c-l) 
2 t 2c + 2 



i 2 A x + 



g 2 c 2 
2t 2 



t~ 2 \i 2 - l)A a 



«9 2 



.9 2c ... d g 2 c (3c - 1) 



tx J -. 

t OX 3 



dx l dx 3 



2t 2 



+ 



£ 2 c (1 — c) . 



<9 2 



<9 2 



2t 2 v y <9t 2 2t 2c + 2V y 2t 2 
We see that even in this simple case, the sectional and curvature diffusion differ tangibly, 
apart from the fact that the range of values of c for which they are defined are different. 



7 References 

[A] Angst J. Etude de diffusions a valeurs dans des varietes lorentziennes. 

Thesis of Strasbourg University, 2009. 
[A-C-T] Arnaudon M., Coulibaly K.A., Thalmaier A. Brownian motion with respect to 
a metric depending on time ; definition, existence and applications to Ricci flow. 
C. R. Acad. Sci. Paris 346, 773-778, 2008. 

[B] Bailleul I. A stochastic approach to relativistic diffusions. 

Ann. I.H.P, vol. 46, n°3, 760-795, 2010. 
[B-E] Beem J.K., Ehrlich P.E. Global Lorentzian Geometry. 

Marcel Dekker, New York and Basel, 1981. 
[De] Debbasch F. A diffusion process in curved space-time. 

J. Math. Phys. 45, n° 7, 2744-2760, 2004. 
[Du] Dudley R.M. Lorentz-invariant Markov processes in relativistic phase space. 

Arkiv for Matematik 6, n° 14, 241-268, 1965. 



32 



[El] Elworthy, D. Geometric aspects of diffusions on manifolds. 

Ecole d'Ete de Probability de Saint-Flour (1985-87), Vol. XV-XVII, 277-425. 
Lecture Notes in Math. Vol. 1362. Springer, Berlin, 1988. 
[Em] Emery, M. On two transfer principles in stochastic differential geometry. 

Seminaire de Probability (1990), Vol. XXIV, 407-441. Springer, Berlin, 1990. 
[F] Franchi J. Relativistic Diffusion in Godel's Universe. 

Comm. Math. Physics, vol. 290, n°2, 523-555, 2009. 
[F-LJ] Franchi J., Le Jan Y. Relativistic Diffusions and Schwarzschild Geometry. 

Comm. Pure Appl. Math., vol. LX, n°2, 187-251, 2007. 
[H-R] Hall G.S., Rendall A.D. Sectional Curvature in General Relativity. 

Gen. Rel. Grav., vol. 19, n°8, 771-789, 1987. 
[H] Harris S.G. A Characterization of Roberts on- Walker Spaces by Null Sectional 

Curvature. Gen. Rel. Grav., vol. 17, n° 5, 493-498, 1985. 
[H-E] Hawking S.W., Ellis G.F.R. The large-scale structure of space-time. 

Cambrige University Press, 1973. 
[Hs] Hsu E.P. Stochastic analysis on manifolds. 

Graduate studies in Mathematics vol. 38, A.M.S., Providence 2002. 
[I-W] Ikeda N., Watanabe S. Stochastic differential equations and diffusion processes. 

North-Holland Kodansha, 1981. 
[K-N] Kobayashi S., Nomizu K. Foundations of Differential Geometry. 

Interscience Publishers, John Wiley & Sons, New York-London-Sydney, 1969. 
[L-L] Landau L., Lifchitz E. Physique theorique, tome II: Theorie des champs. 

Editions MIR de Moscou, 1970. 
[M] Malliavin, P. Geometrie differentielle stochastique. 

Les Presses de l'Universite de Montreal. Montreal, 1978. 

A. M.S. Classification : Primary 58J65, secondary 53C50, 83C10, 60J65. 

Key Words : Covariant relativistic diffusions, General relativity, Lorentz manifold, Rie- 
mann curvature tensor, Stochastic flow, Robert son- Walker manifolds. 

• Jacques FRANCHI : jacques.franchi@math.unistra.fr 

IRMA, Universite de Strasbourg et CNRS, 7 rue Rene Descartes, 67084 Strasbourg cedex, France ; 
beneficiaire d'une aide de l'Agence Nationale de la Recherche, n° ANR-09-BLAN-0364-01. 



33 



• Yves LE JAN : yves.lejan@math.u-psud.fr 

Universite Paris Sud 11, Mathematiques, Batiment 425, 91405 Orsay, France. 



34 



